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^. , In this paper we study the Lyapunov stability and Hopf bifurcation in a 

system coupling a Watt-centrifugal-governor with a steam-engine. Sufficient 
conditions for the stability of the equilibrium state in terms of the physical 
parameters and of the bifurcating periodic orbit at most critical parameters 
on the bifurcation surface are given. 
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^ 1 Introduction 

The Watt centrifugal governor is a device that automatically controls the 
speed of an engine. Dating to 1788, it can be taken as the starting point 
for automatic control theory (see MacFarlane [5. and references therein). In 
this paper the system coupling the Watt-centrifugal-governor and the steam- 
engine will be called simply the Watt Governor System (WGS). 

Landmarks for the study of the local stability of the WGS are the works 
of Maxwell jB] and Vyshnegradskii jHj . A simplified version of the WGS local 
stability based on the work of Vyshnegradskii is presented by Pontryagin [JJ. 
A local stability study generalized to a more general Watt governor design 
was carried out by Denny |2j. 

Enlightening historical comments about the Watt governor local mathe- 
matical stability and oscillatory analysis can be found in MacFarlane [5] and 
Denny [2]. There, as well as in [JJ, we learn that toward the mid XIX cen- 
tury, improvements in the engineering design led to less reliable operations in 



the WGS, leading to fluctuations and oscillations instead of the ideal stable 
constant speed output requirement. The first mathematical analysis of the 
stability conditions and subsequent indication of the modification in design 
to avoid the problem was due to Maxwell [U] and, in a user friendly style 
likely to be better understood by engineers, by Vyshnegradskii [S|. 

From the mathematical point of view, the oscillatory, small amplitude, 
behavior in the WGS can be associated to a periodic orbit that appears from 
a Hopf bifurcation. This was established by Hassard et al. in >3J and Al- 
Humadi and Kazarinoff in pQ. Another procedure, based in the method of 
harmonic balance, has been suggested by Denny [2] to detect large amplitude 
oscillations. 

In this paper we present a simple way to understand the Hopf bifurca- 
tion in a WGS, which is more general than that presented by Pontryagin [7], 
Al-Humadi and Kazarinoff [1J and Denny j2j. We believe that our approach 
has the advantage of connecting the physical parameters of the system to 
the stability of the stationary equilibrium point and the bifurcating periodic 
orbit. It also permits a neater geometric synthesis of the bifurcation analysis 
based on the algebraic expression and geometric location of the curve — the 
codimension 2 Hopf points — characterizing the stability versus the instabil- 
ity at the critical weak focal equilibria on the critical surface of parameter 
values — the codimension 1 Hopf points — leading to the bifurcation of pe- 
riodic orbits, whose stability depend on the side of the curve at which the 
parameters cross the surface. 

This paper is organized as follows. In Section El we introduce the differ- 
ential equations that model a general WGS . The stability of the equilibrium 
point of this model is analyzed and a general version of the stability condition 
is obtained and presented in the terminology of Vyshnegradskii. The codi- 
mension 1 Hopf bifurcation for the WGS differential equations is studied in 
Section El An expression — neater than that found in the current literature 
- which determines the sign of the first Lyapunov coefficient is obtained. 



Sufficient conditions for the stability for the bifurcating periodic orbit are 
given. Concluding comments are presented in Section 0] 

2 A general Watt governor system 

2.1 General differential equations 

The WGS studied in this paper is shown in Fig. [TJ There, tp G (0, |) is 
the angle of deviation of the arms of the governor from its vertical axis S\, 
Q G [0, oo) is the angular velocity of the rotation of the engine flywheel D, 9 
is the angular velocity of the rotation of S\, I is the length of the arms, m is 
the mass of each ball, if is a special sleeve, T is a set of transmission gears 
and V is the valve that determines the supply of steam to the engine. 

The differential equations of our model which generalize those found in 
Pontryagin [7j, p. 217, are given by 
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(2) 
where r is the time, ip — -£, g is the standard acceleration of gravity, 
6 = s(Q), s is a smooth non-negative increasing function with s(0) = 0, 
called the transmission function, h is a non-negative increasing function with 
h(0) = that represents the frictional force of the system, I is the moment 
of inertia of the flywheel D and M is a smooth decreasing function of the 
angle <p. This function depends upon the design of the WGS and determines 
the supply of steam to the engine through the valve V. It measures the effect 




Figure 1: General- Watt-centrifugal-governor-steam-engine system. 

that the dynamics of the governor has on the engine. Denny in |5J has also 
considered this function for the torque. The torque F due to the load is given 
by F = M(</?*), where (p* is an equilibrium angle. 

The standard Watt governor differential equations in Pontryagin [7j, p. 
217, 
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are obtained from (£[]) taking 



s(Q) —eft, h(ip) = b tp, M(p) = /i cos <p, c, b, \i > 0. 



(4) 



2.2 Stability analysis of the equilibrium point 



Let P = (foj 0, fi ) be an equilibrium point of ((TJ). The Jacobian matrix of 
/ at Pq has the form 
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For the sake of completeness we state the following lemma whose proof 
can be found in j2], p. 58. 

Lemma 2.1 The polynomial L(X) = Po\ 3 +Pi\ 2 +P2X+P3, Po > 0, with real 
coefficients has all roots with negative real parts if and only if the numbers 
Pi,P2,P3 are positive and the inequality pip 2 > P0P3 is satisfied. 



Theorem 2.2 If 



h'(0) > 



2 m s'(Qq 



M'((p ) cottpo 



(6) 



I s(Q j 

then the WGS differential equations J7J) have an asymptotically stable equi- 
librium point at Pq. If 

, .. „ 2 m s'(Qn) , , 
< ti{0) < — -y-^ M'{(p ) cot^o 



I s(fio) 



then P n is unstable. 



Proof. The characteristic polynomial of Df (Pq) is given by p(\), where 



(X) =- A 3 + h '^ X 2 + 9 Sln2 ^° x 2g M '^ S '^ sin ^° 



A 
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As h'(0) > 0, s'(Q ) > 0, s(Q ) > and M'{<p ) < the coefficients of -p(X) 
are positive. Thus a necessary and sufficient condition for the asymptotic 



stability of the equilibrium point Pq, as provided by the condition for one 
real negative root and a pair of complex conjugate roots with negative real 
part, is given by (JBJ), according to Lemma 12. 11 



In terms of the WGS physical parameters, condition (0) is equivalent to 

h'{0) I 



m 

where 

7/ 



V >1, (8) 
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s(fto) 

tany3o 



2 M'(<po) s'(ft ) 

is the non-uniformity of the performance of the engine which quantifies the 
change in the engine speed with respect to the load (see |7j, p. 219, for more 
details). The rules formulated by Vyshnegradskii to enhance the stability 
of the WGS follow directly from (jHJ). In particular, the interpretation of (0) 
is that a sufficient amount of damping h'(0) must be present relative to the 
other physical parameters for the system to be stable at the desired operating 
speed. The general condition (JHJ) is equivalent to the original conditions given 
by Vyshnegradskii (see [7j, p. 219). 

In next section we study the stability of Pq under the condition 

, ,. . 2 m s'(Qn) , ,.. s 

h'(0) = r 4?7T M ^ cot ^> ( Q ) 

I s(il ) 

that is, on the surface — the Hopf surface — complementary to the range of 
validity of Theorem 12.21 

3 Hopf bifurcation analysis 

For the analysis carried out here we take (0) with s(Q) = c Q, where c > 
is a constant transmission ratio and h(ip) = bip, where b > is a constant of 
the frictional force of the system. 
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After the following change in the coordinates and the time 

x = (p } y = J-4>, z = cJ-Q, r=J~t, (10) 

the differential equations (JIJ can be written as 



:n) 



or equivalently as x' = /(x) where 

x= (x,y,z) e(o^)xRx [0,oo), 

/(x) = (y, z 2 sin x cos a; — sinx — e y,T(x) — G) , (12) 

b I r^, x c I ,. . _ cl F 
e = -J-, T(x) = — Ms, G = ——. 
m y g Ig Ig 

Here e is considered to be a changing parameter; so that the differential equa- 
tions (jll)) (or in its equivalent vectorial form (fT2|) ) can in fact be regarded 
as a family of one-parameter families of differential equations, dependent on 
the functional parameter T. 

In this section we will analyze the stability at 



P = (x ,y ,zo) = ( arccos/3,0, J- J 



(13) 



under the condition Q, which now writes as 



UJ 

where (3 = cos x and 



2 V(x ), (14) 



W0 = \l n ■ ( 15 ) 



3.1 Generalities on Hopf bifurcations 

The study outlined below is based on the approach found in the book of 
Kuznetsov |3j, pp 177-181. 

Consider the differential equations 

x' = /(x,/i), (16) 

where x e M 3 and \x e W 1 is a vector of control parameters. Suppose (|16|) 
has an equilibrium point x = x at /i = /i and represent 

F(x) = /(x,/i ) (17) 

as 

F(x) = Ax+ ~ B(x,x) + ~ C(x,x,x) + 0(||x|| 4 ), (18) 

2 o 

where A = / x (0,/xo) an d 
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C<(x,y,z)= V — ' XjykZi, (20) 

J^i <^ % % e=o 

for i = 1,2, 3. Here the variable x — xo is also denoted by x. 

Suppose (x , Ho) is an equilibrium point of (fTo)) where the Jacobian matrix 
^4 has a pair of purely imaginary eigenvalues A 2i3 = ±iou , uj q > 0, and no 
other critical (i.e., on the imaginary axis) eigenvalues. 

The two dimensional center manifold can be parametrized by w G K 2 = 
C, by means of x = H(w, w), which is written as 

H(w,w) = wq + wq+ > — — ■ h jk w^w k + 0(|u>| 4 ), 

2<j+k<3 J 



with hjk € C 3 , hjk = h 



hi 



Substituting these expressions into (|16|) and (|T5jl we have 

H w (w, w)w' + H a (w, w)w' = F(H(w, w)). (21) 

Let p, q G C 3 be vectors such that 

3 

Aq = iu q, A T p = -iu p, (p, q) = ^pi qi = 1. (22) 

i=X 

The complex vectors h^ are to be determined so that equation (|21|) writes 
as follows 

w' = IUJqW H — G2iu>|w| 2 + 0(|w| 4 ), 

with G 2 i e C. 

Solving the linear system obtained by expanding (|2*T|) . the coefficients of 
the quadratic terms of (|T7jl lead to 

/H^-A^fi^g), (23) 

/i 20 = (2^ / 3 -A)- 1 J B(g,g), (24) 

where I3 is the unit 3x3 matrix. 

The coefficients of the cubic terms are also uniquely calculated, except 
for the term w 2 w, whose coefficient satisfies a singular system for /i 2 i 

(iuoh ~ A)h 21 = C{q, q, q) + B{q, h 20 ) + 2B(q, h lt ) - G 21 q, (25) 

which has a solution if and only if 

(P, C(q, q, q) + B(q, h 20 ) + 2B(q, h u ) - G 21 q) = 0. 

Therefore 

G21 = (p,C(q 1 q 1 q)+B(q,(2tuJoh-A)- 1 B(q,q))-2B(q 1 A~ 1 B(q,q))}, (26) 

and the first Lyapunov coefficient l\ - which decides by the analysis of third 
order terms at the equilibrium its stability, if negative, or instability, if pos- 
itive - is defined by 

l 1 = J-ReG 21 . (27) 

Z UJq 
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A H op f point (x ,/io) is an equilibrium point of (J16JI where the Jacobian 
matrix A has a pair of purely imaginary eigenvalues A2,3 = iicuo, ^o > 0, 
and no other critical eigenvalues. At a Hopf point, a two dimensional center 
manifold is well-defined, which is invariant under the flow generated by (|16|) 
and can be smoothly continued to nearby parameter values. 

A Hopf point is called transversal if the curves of complex eigenvalues 
cross the imaginary axis with non-zero derivative. 

In a neighborhood of a transversal Hopf point with 1^0 the dynamic 
behavior of the system ()16j). reduced to the family of parameter-dependent 
continuations of the center manifold, is orbitally topologically equivalent to 
the complex normal form 



w' — (7 + iuj)w + liw\w\ 



(28) 



w G C, 7, uj and l\ are smooth continuations of 0, uoq and the first Lyapunov 
coefficient at the Hopf point |3]. When l\ < (/1 > 0) a family of stable 
(unstable) periodic orbits appears found on this family of center manifolds, 
shrinking to the equilibrium point at the Hopf point. 

From ()12|) write the Taylor's expansion (fTHj) of /(x). Define a\ = T'(xq), 
a 2 = T"(x ) and a 3 = T'"(x ). Thus 

/ 1 \ 



A 



-ul -e c 2(3uj 



(29) 



\ ai / 

and, with the notation in (|T8|) we have 

F(x) - Ax = (0,F 2 (x) + 0(||x|| 4 ), ^ x 2 + I x 3 + 0(||x|| 4 ) ) , (30) 
where 



FoU) 



2 ° 



(3x 2 + u; ^ 2 z 2 + 



2(2(3 2 - 1) 



x z + 



4-7/3 2 



x 



VP 6/3 

AujqPx 2 z + (2(3 2 -l)xz 2 . 
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From (|29j) the eigenvalues of A are 

Ai = -e c , \ 2 = iuj 0j X 3 = -iuJo. (31) 

The eigenvectors q and p satisfying (}2"2"j) are respectively 



and 

The main result of this section can be formulated now. 

Theorem 3.1 Consider the family of differential equations Ml)) - The first 
Lyapunov coefficient at the point HH\) for parameter values satisfying ([?7P is 
given by 

2 ((1 - /3 2 ) 4 + 5/? 4 (l - f3 2 ) 2 aj + 4/5X) 
where 



h{f3,a u a 2l a 3 ,E c ) = - :777: ^ on4 : ^77; 2 — — ^ , (34) 



H(/3, 01, oa, a 3 ) = /5 2 2/3 7 a? + /3 3 a^2 - 3/3 2 + 5/3 4 ) + 4/3(1 - (3 2 fa 2 2 + 
/3(1 - /? 2 ) (a 2 (9(1 - /? 2 ) 2 + 2/? 4 a 2 ) + 2a 3 ( - /3 4 a? - (1 - /? 2 ) 2 a x ) ) + 
a 2 V / ^ 3 ^( - /? 4 af(l + 5/3 2 ) - oi(l - /? 2 ) 2 (5 + 3/3 2 )) 



Proof. The proof depends on preliminary calculations presented below. 
From (HHJ), (HH, flUl) and ((30j) one has 

5(x,y) = (0,S 2 (x,y),a 2 a:il/i), (35) 

where 

2 (2 (3 2 — I) 
5 2 (x, y) = -3w 0V ^i 1 ?/ 1 + 2 wo /9 3/2 ^3 2/3 H 7=5 (a?i J/a + ^3 Vi) , 
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C(x, y, z) = (0, C 2 (x, y, z), a 3 x x y x z x ) , 

where 

4-7/3 2 
C 2 (x, y, z) = — xi t/i ^i - 8^0 P {xi V\ z 3 + x x 2/3 z x + x 3 y x z x ) + 



(36) 







2 {2j3 2 - 1) (xi y 3 z 3 + x 3 y x z 3 + x 3 y 3 z x ) . 



Referring to the notation in (USD, (J2EJ), (JHD, (ED and (J22D one has 



C(?,S,s)= (o, 

where 
B 2 (q,h u ) 



A(3 2 a x 



(37) 

(38) 

-z o 3 ) , (39) 

(40) 

(6/? - £ 2 ) (2 2 - 4z/? 2 + /fe^oK - 2 V / /?( 2i(l + (3 2 )uj q + 

e c (-2 + /3(4/3 + ^ 2 ))')a 2 



2v^ 
g/3 2 gc a;o + i {(5{U(5 2 - 8) - e 2 (2/3 2 - 1)) 

B(q,h lx ) = [0,B 2 {q,h lx ),i^ 



and 

where 

B 2 (q,h 20 



B{q, h 20 ) = (0, £ 2 (g, h 20 ), B 3 (q, h 20 )) 



(41) 



/? 12^ 2 (2/3 2 -lVo + 



4/5 3 cj ((2e c + 3z^o)cuo + /5a 1 ) 
6i/?(2/? 2 - l)e 2 w - fi 2 e\ul + e c (8- 32(3 2 + 32/3 4 - 6/? 3 cj 2 ) ) oi + 



2u 1 - 18i/? 2 (2/3 2 - l)e c cj - i/3(2/? 2 - l)e°w + 18/Tu; 2 + 
(-4 + 16/3 2 - 16/3 4 + 3/3 3 w 2 ) + /3 5/2 (6i(3/3 2 - l)u + 2i(3u e 2 c + 



3e c (-2 + 4p 2 -pu 2 ))a 2 



12 



and 

a 2 (/3u (6/3 + el) + i ((4 - 8/3 2 )e c + 2/3 5 / 2 a 2 )) 
3lg ' 20j 2/3 3 /2(2 £fM) + /3a 1 + 3^ 2 ) 

The first Lyapunov coefficient is given by (|27|) . From ()33|) and (|39|) one 



has 



Re(p,C(g,g,g)>- 4/3 2 K 2 + . 2 ) ' (42) 



From (J33J) and (|3UJl one has 

1 



Re<p, 2fl(g, /m)) 



e c (6/3 - £ 2 )(-2 + 4/3 2 + /3wj})ai + 

(43) 



2y/pu a 2 (e c (4 - 2/3 2 - /3u; 2 ) + 4/3 5 / 2 a 2 ) 



From ()33|) and (|41jl one has 



R / R(~ h \\ = $(p,ai,a 2 ,e c ,uj ) 

{P ' W ' 20>) 8/33 (w g + £ 2 } ( 4 £cW 2 + ^4 + p ai{4£cUo + ^0) ' l ] 

where 

0(/3, en, a 2 , e c , w ) = 4(2/3 2 - V)e\ (10 + /3(e 2 - 26/3)) a; 2 + /3 U f 2(1 - /? 2 ) 

e 2 (/3(42/3 + e 2 ) - 24H - /3(36/3(3/3 2 - 1) + 6(1 + 4/3 2 )e 2 + 5/3^)^ - 

/?e c (2(2/3 2 - l)(/3(2/3 + e 2 ) - 4) + /? 2 (6/3 + £ c 2 Vo>i) + 2/? 3/2 ^o 

(18(1 - 5/3 2 )^ + 4^(38 - 5/3(8/3 + 3u$) + /3 £c (4 + 10/3 2 - 3/3wjj)ai) a 2 + 

8/3V)(5£ c w + /3ax)a 2 ) - 6/3e^ (6/3(1 + 3/3 2 ) + (7/3 2 - l)e 2 ) . 



Substituting flEJ) and (JUJ) into (g2), (02)) and (jSJ) and the results into (J2SD 
and (|2*7|). the theorem is proved. 
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Proposition 3.2 Consider the family of differential equations ill)) regarded 
as dependent on the parameter e. The real part, 7, of the pair of complex 
eigenvalues verifies 



l'{e c ) 



1 wk 



<0. 



2ul + el 
Therefore, the transversality condition holds at the Hopf point. 



(45) 



Proof. Let A(e) = A2,3(e) = 7(e) ± iuj{e) be eigenvalues of A[e) such that 
7(e c ) = and u(e c ) = u , according to (|3*T|). Taking the inner product of p 
with the derivative of 

A(e)q(e) = X(e)q(e) 



at e = e c one has 






/ dA 

\ ' de 


q\=i(e c )±u\e c ). 

£=£c 1 


Thus the transversality condition is given by 


.. . 1 dA 
7 (e c ) = Re {p, -^ 


£=e c / 


As 

de 


q = (0, - 


^o,0), 



the proposition follows from a simple calculation. 



(46) 



Theorem 3.3 Define Sp = {(02,03)! a 3 + K 2 (f3)a 2 > 0}, where 

If (02,03) G S'/j then the one-parameter family of differential equations Ml\) 
has a Hopf point at Pq for e = e c . Furthermore this Hopf point at Pq is 
asymptotically stable and for each e < e c , but close to e c , there exists a stable 
periodic orbit near the unstable equilibrium point Pq. See Fig. [H 
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Proof. Since the denominator of (|34|) is positive the sign of the first Lya- 
punov coefficient is determined by the sign of R{{3, Oi, a 2 , a 3 ). Rewrite this 
expression in the following way 
R((3, a x , a 2 , a 3 ) 



/?(l-/3 2 )(a 2 (9(l-/? 2 ) 2 + 2/3V 2 ) 



2(3 7 al + /3 3 a 4 (2 - 3/3 2 + 5/3 4 ) + 4/5(1 - {3 2 fa\ + 

-2/3 5 a?(l-/3 2 )(a 3 + ^ 1 (/5)a 2 )- 
2a 1 /3(l-/5 2 ) 3 (a 3 + ^ 2 (/3)a 2 ), 



where 



K.iP) = ^±ML > 0, 
2(3^T 2 



and K 2 (f3) is as above. Since 2 — 3/3 2 + 5/3 4 is positive, the sum into the 
bracket is positive. A simple calculation shows that 

K 2 (/3) > KM, 

for all (3 G (0, 1). If (a 2 , a 3 ) G Sp then 

-2/3 5 a 3 (l - /3 2 )(a 3 + K^a,) - 2(1,(3(1 - (3 2 ) 3 (a 3 + K 2 {f3)a 2 ) > 0, 

since 01 < 0. This implies that R{(3, a%, a 2 , o 3 ) > and therefore the first 
Lyapunov coefficient is negative. See Fig. |21 

■ 

Remark 3.4 27ie family ill)) has been studied in Denny J$j, focusing the 
Lyapunov stability of the equilibrium point ilfy) . outside the bifurcation sur- 
face (fl^l ). The stability conditions in [7|/ and J$ have been extended in The- 
orem \2. 6 2i to a more general WGS model given in (QJj. For this more general 
model, however, the study of the codimension one Hopf bifurcation involves 
much longer calculations. For this reason, the analysis performed here has 
been restricted to the case proposed by J$/. See Theorem \3.§. The calculations 
for the general case in (QJ) ; while too long to put in print, can be handled by 
Computer Algebra. The special case of the equations of Pontryagin ^ will 
be treated in next subsection. 
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a 3 . 




\ a 2 



Figure 2: The region Sp . 

3.2 Hopf Bifurcations in Pontryagin equations 

Below we specialize the calculations above to the case of the Eq. (J3J) for the 
WGS in the form presented by Pontryagin [7] also treated in p. 

With the change in the coordinates and time fflUj) the standard Watt gov- 
ernor differential equations (J3j) are obtained from (jlljl with M(x) = ficosx 
and have the form 



x = y 

y' = z 2 sin x cos x — sin x — e y 

z' — a (cosx — j3) 



(47) 



where a > 0, (3 G (0, 1) and e > are given by 




J 1 p 

9 1 V 



(48) 



The differential equations (J47J1 have an equilibrium point located at Pq 
given by Eq. (|T3j) . 
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A necessary and sufficient condition for the hyperbolic asymptotic stabil- 
ity of the equilibrium point Pq is e > 2af3 3 ^ 2 , according to Theorem 12 .21 For 
< £ < 2 a /3 3//2 , the equilibrium point Pq is unstable. Now we analyze the 
stability of P as e c = e(/3, a) = 2 a f3 3 ^ 2 obtained from (|T3j) . 

Theorem 3.5 Consider the family of differential equations fl^Tj ). The first 
Lyapunov coefficient on the Hopf surface e = e c = 2 a /3 3//2 is given by 



j (B , af3 2 ^T^(3 + (a 2 - 5)/? 2 + a 4 /? 6 ) 

llP ' a ' £cJ 2 (1 - (3 2 + a 2 /3 4 ) (1 - (3 2 + Aa 2 ^) 



(49) 



Proof. The proof follows from ([34J1 with 

oi = -ayl — (3 2 , a 2 = —a (3, a 3 = ayl — /3 2 . 



Theorem 3.6 // 

(/3, a ) = 3 + (a 2 - 5)(3 2 + a 4 /? 6 (50) 

is different from zero then the family of differential equations J^7[ ) /ias a #op/ 
point at Pq for e c = 2a /3 3//2 . 



Proof. From ()45|) the transversality condition is satisfied. Therefore a suf- 
ficient condition for being a Hopf point is that the first Lyapunov coefficient 
li(j3,a,e c ) 7^ 0. But from (J4T?|) it is equivalent to g((3,a) ^ 0. The theorem 
is proved. 



The following theorem summarizes the results of this subsection. 

Theorem 3.7 If (/?, a, e c ) G S U U then the family of differential equations 
(OTp /icts a Hopf point at P . If (/3, a, e c ) G S 1 t/ien t/ie Hopf point at P is 
asymptotically stable and for each e < e c , but close to e c , there exists a stable 
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periodic orbit near the unstable equilibrium point Pq. If (/3,a,e c ) G U then 
the Hopf point at Pq is unstable and for each e > e c , but close to e c , there 
exists an unstable periodic orbit near the asymptotically stable equilibrium 
point Pq. 

Corollary 3.8 Consider the family of differential equations fl^7| ). If a > 1 
then the equilibrium point Pq is asymptotically stable for e = e c and for all 
< /3 < 1. Therefore for each e < e c , but close to e c , there exists a stable 
periodic orbit near the unstable equilibrium point Pq. 

Proof. The proof is immediate from Theorem 13.71 In fact, if a > 1 then 
(P, a, e c ) G S for all < (3 < 1. Thus the first Lyapunov coefficient h(@, a, e c ) 
is negative. 
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Figure 3: Signs of the first Lyapunov coefficient. 



Remark 3.9 Equation \5U\) gives a simple expression to determine the sign 
of the first Lyapunov coefficient \49\) . Its graph is illustrated in Fig. where 



the signs of the first Lyapunov coefficient are also represented. This gives an 
analytic corroboration to the observation in J3j/, p. 155, made on the basis of 
numerical testing: 

"For a > 1, our computations indicate that li((3,a,e c ) < for all values 
of (3, 0<(3<1." 

This matter has been addressed in £7J/ ; providing a longer expression to 
decide the sign first Lyapunov coefficient than the one given in fijjfy . 

The curve l\ = divides the surface of critical parameters e c = 2 a /3 3//2 
into two connected components denoted by S and U where l\ < and l\ > 
respectively. See Fig. 

3.3 Examples of stability and instability in the WGS 

The following examples give a piece of information about the stability of the 
bifurcating periodic orbit of the family ()47|) . First in terms of the angle of 
deviation (p. 

Example 3.10 If the equilibrium angle of deviation 

ifio > 39.23° (51) 

then the equilibrium point Pq is asymptotically stable for e = e c and for all 
a > 0. Therefore for each e < e c , but close to e c , there exists a stable periodic 
orbit near the unstable equilibrium point Pq. 

In fact, if <p > 39.23° then (3 = cos^ < 0.7746. Therefore l x {e c ) < for all 
a > 0, as follows from Fig. El 

Now in terms of the normalized speed. 

Example 3.11 If the normalized engine speed zq satisfies 

z > 1.1362 (52) 

then the equilibrium point Pq is asymptotically stable for e = e c . Therefore 
for each e < e c , but close to e c , there exists a stable periodic orbit near the 
unstable equilibrium point P . 
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In fact, if z Q > 1.1362 then ip Q > 39.23° since cos<^ = V^o- The proof 
follows from the Example 13.101 

In terms of the normalized parameter (3 in (|48j). we have: 

Example 3.12 If 

0.7746 < < 1 

and 



v/20/3 4 - 12/3 2 + 1-1 
< a < = 

V2P 2 

then the equilibrium point Pq is unstable for e = e c . Therefore for each e > e c , 
but close to e c , there exists an unstable periodic orbit near the asymptotically 
stable equilibrium point Pq. 

This follows from Fig. El In fact, the level curve g(/3,a) = 0, where g(j3,a) 
is defined in (150(1 , is the graph of the function 



J20(3 4 - 12/3 2 + 1-1 
h{(3) = ^-L . 

4 Concluding comments 

The historical relevance of the Watt governor study as well as its importance 
for present day theoretical and technological control developments - going 
from steam to diesel and gasoline engines - have been widely discussed by 
MacFarlane jH] , Denny [2] and Wellstead - Readman [2] among others. 

In this paper the original stability analysis due to Maxwell and Vyshne- 
gradskii of the Watt Centrifugal Governor System -WGS- has been revisited 
and extended to the following situations: 

1. A more general torque function of the engine has been considered. This 
function determines the supply of steam to the engine via the valve V 
and is related to the design of the WGS. It measures the effect of the 
governor on the engine; 
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2. A more general transmission function has been studied. This function 
determines the relation between the angular velocities of the engine and 
the governor axis. It measures the effect of the engine on the governor; 

3. A more general frictional force of the system has been modelled. 

In Theorem 12.21 we have extended the stability results presented in Pon- 
tryagin [7j and Denny j2|. 

Concerning the bifurcations of the WGS this paper deals with the codi- 
mension one Hopf bifurcations in the Watt governor differential equations. 
The main results are: 

1. In Theorem 13.31 we give sufficient conditions for the stability of the 
periodic orbit that appears from the Hopf point for the general Watt 
governor differential equations (jllj) ; 

2. In Theorem 13 .71 we have extended and provided a neat geometric syn- 
thesis of the Hopf stability analysis performed by Hassard et al. jH] 
and Al-Humadi and Kazarinoff Ij. This identifies the codimension one 
Hopf points and the location of its complement codimension 2 Hopf 
points along a curve in Fig. El 
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